Universal Spin Response in Copper Oxide Materials 
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The spin response in the copper oxide materials at finite temperatures in the underdoped and 
optimal doped regimes is studied within the framework of the fermion-spin theory. The integrated 
dynamical spin structure factor is almost temperature independent, the integrated susceptibility 
shows the particularly universal behavior, and the spin-lattice relaxation time is weakly temperature 
dependent, which are consistent with experiments and numerical simulations. 
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After ten years of intense experimental and theoreti- 
cal studies of the copper oxide superconductors, there is 
now a consensus that these materials should be described 
as strongly correlated electron systems 1 , since all of the 
copper oxide materials have in common the existence of 
a perovskite parent compound which is insulating and 
has the antiferromagnetic long-range-order (AFLRO), 
and changing the carrier concentration by ionic sub- 
stitution or increase of the oxygen content turns these 
compounds into correlated metals leaving short range 
antiferromagnetic correlations still intact 1 ' 2 . The short 
range antiferromagnetic correlations results in several pe- 
culiar physical properties of the copper oxide materi- 
als: they are responsible for the nuclear magnetic reso- 
nance (NMR) and nuclear quadrupole resonance (NQR) , 
and especially for the temperature dependence of the 
spin-lattice relaxation rate 1 ' 3 . A series of the neutron- 
scattering measurements 4-6 on the copper oxide materi- 
als La2- a; Sr a ;Cu04 and YB^CuaOe+z show that there 
is an anomalous temperature T dependence of the spin 
fluctuations near the antiferromagnetic zone center in the 
underdoped and optimal doped regimes, and the low- 
frequency dynamical susceptibility in the optimal doped 
regime follows a surprisingly simple scaling function as 
X"M arctan(w/T). The NMR and NQR spin-lattice 
relaxation time T± is weakly T dependent. These un- 
usual magnetic properties of the copper oxide materials 
suggest that the normal-state can not be described by 
the conventional Fermi-liquid theory. 

As emphasized by many researchers 7 , the essential 
physics of the copper oxide materials is contained in 
the doped antiferromagnet, which may be effectively de- 
scribed by the two-dimensional (2D) t-J model acting on 
the space with no doubly occupied sites. The t-J model 
is reduced as the Heisenberg model in the undoped case. 
In spite of its simple form the t-J model proved to be 
very difficult to analyze, analytically as well as numer- 
ically, because of the electron single occupancy on-site 
local constraint. The local nature of the constraint is of 
prime important, and its violation may lead to some un- 
physical results 8 . Recently a fermion-spin theory based 
on the charge-spin separation is proposed 9 ' 10 to incorpo- 



rated this constraint, where the electron on-site local con- 
straint for single occupancy is satisfied even in the mean- 
field approximation (MFA). The magnetic in the un- 
doped parent compounds is now quite well understood 1 : 
here, the system of interacting localized Cu 2+ spins is 
well described by the 2D Heisenberg model, and then 
it is clearly of great interest to investigate in detail the 
crossover from the rather conventional local moment sys- 
tem at zero doping to the electronic state that forms the 
basis for the high-temperature superconductivity. There- 
fore in this paper, we only study the spin dynamics of the 
copper oxide materials within the fermion-spin theory in 
the underdoped and optimal doped regimes. According 
to the fermion-spin formulism 9 ' 10 , the electron operators 
can be decomposed as d\ = h\S~ and Cjj = h\Sf , with 
the spinless fermion operator hi keeps track of the charge 
(holon), while the pseudospin operator 5, keeps track of 
the spin (spinon). Within the fermion-spin theory, it 
has been shown 11 that AFLRO vanishes around doping 
6 = 5% for the reasonable value of the parameter t/J = 5. 
The mean-field theory in the underdoped and optimal 
doped regimes without AFLRO has been developed 10 , 
where the mean-field order parameters are defined as 
X = (S i S i+J1 ) — (S i S i+rj ), Xz = (S*S* +r) ), c = 

{^/ z2 )H, h r,<{si +11 s l+1 y), c z = (^/z 2 )Y /r ,^(Si +v sf +71 ,), 

and (f> — {h\hi +ri ) with f) = ±i, ±y, and Z is the number 
of nearest neighbor sites. In this case, the low-energy 
behavior can be described 10 by the effective Hamiltonian 
H = H t + Hj with 

H t = -tJ2h t hl +v (S+S- +v + Sr S + +v ) + h.c. 

if] 

i 

H.j = Jeff y^\n{ S t Sj+n + S i S t+7 1 ) + S i S i+n}> ( lb ) 
irj 

where J e // = J[(l — <5) 2 — cj) 2 ], and /i is the chemical 
potential which enforce (hi hi) = S. 

In the framework of the charge-spin separation, the ba- 
sic low-energy excitations are holons and spinons. The 
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charge dynamics can be discussed based on the Ioffe- 
Larkin combination rule 12 , however, since the spin fluc- 
tuations couple only to spinons, and therefore no compo- 
sition law is required in discussing the spin dynamics 12 , 
but the strongly correlation between holons and spinons 
still is considered through the holon's order parameters <f> 
entering in the spinon propagator, which means that the 
spinon moves in the background of holons, and the cloud 
of distorted holon background is to follow spinons, there- 
fore the dressing of the spinon by holon excitations is the 
key ingredient in the explanation of the spin dynamics. 
According to Ioffe-Larkin combination rule 12 , we 13 have 
discussed the optical conductivity, Drude weight, and re- 
sistivity of the copper oxide materials in the undcrdopcd 
and optimal doped regimes by considering fluctuations 
around the mean-field solution, where the dominant dy- 
namical effect is due to the strongly spinon-holon inter- 
action in Hamiltonian (1). We believe that this strongly 
spinon-holon interaction also will dominate the spin dy- 
namics within the same regimes. The mean-field spinon 
Green's functions D(°\k, lu) and D^\k,u)) and mean- 
field holon Green's function g^°\k,w) have been given 
in Ref. 10 . In this paper, we limit the holon part to the 
first-order (mean-field level) since some physical prop- 
erties can be well described at this level 10 , and spin 
fluctuations couple only to spinons as mentioned above. 
However, the second-order correction for the spinon is 
necessary for the discussion of the spin dynamics. The 
second-order spinon self-energy diagram from the holon 
pair bubble is shown in Fig. 1. Since the spinon opera- 
tors obey the Pauli algebra, we map the spinon operator 
into the spinless-fermion representation in terms of the 
2D Jordan- Wigner transformation 14 for the formal many 
particle perturbation expansion. After then the spinon 
Green's function in the spinon self-energy diagram shown 
in Fig. 1 is replaced by the mean-field spin Green's func- 
tion D^°\k,u). In this case, we obtain the second-order 
spinon self-energy as, 

£< 2 )(fc, W ) = -(Ztf^Yt^-p+lp'+p+kfBk+r, 
pp' 

x / F 1 (k,p,p') 

\lo + £ p+p / - i p + Uk+p' + i0 + 

F2(k,p,p') \ 

<*> + £p+p' ~€p- Uk+p' + z0+ / ' 

where F 1 (k,p,p') = ra F (£ p+p /)[l - n F (€ P )] + [1 + 
n B (^k+p')}[nF^ P ) - n F (£ p+p /)], 

Fi (k,P, p') = n F (£p+ P > ) [1 - n F (€ P )] - n B (oJk+ P ' ) [n F (£ P ) - 
MW)L 7fc = (l/^)E,e ifc -*», e = l + 2^/J e// , 
B k = ZJ eff [(2eXz + x)lk - (ex + 2 Xz )]/^k, n F (£ k ) 
and nfl(wjt) are the Fermi and Bose distribution func- 
tions, respectively, the mean-field holon excitation spec- 
trum £fc = 2Z\tj k + an d the mean-field spinon ex- 
citation spectrum uj k is given in Ref. 10 . Then the full 
spinon Green's function is obtained as D^ 1 (k,uj) — 
£)(°)- 1 (fc ; w )_ Y,f\k, uj). Since the local constraint of the 



t-J model has been treated exactly in the previous mean- 
field theory 10 , and it is natural satisfied in the above per- 
turbation expansion based on this mean-field theory. 

We are now ready to discuss the spin dynam- 
ics. The dynamical spin response, as manifested 
by the dynamical spin structure factor S(k,u)) and 
the susceptibility x(fc,w), are given as S(k,oj) = 
Rej^dte lut (S+{t)S^{Q)) = 2ImD(fc, w)/(l - e~^ u ) and 
X"(k,u) = (1 - e-^)S(k,uj) = 2ImD(fc, w). The prop- 
erties of S(k,u) and x"{k,u) in different k directions 
have been discussed 15 , and the results showed that there 
is the anomalous temperature T dependence of the spin 
fluctuations near the antiferromagnetic point Q — (tt, it). 
In this paper we are interested in the universal behav- 
ior of the integrated dynamical response. The integrated 
dynamical spin structure factor and integrated suscepti- 
bility are expressed as, 

S(u>) = S l (cj) + S l (-cj) = (1 + e~ pu )S L {u), 

k 

and 

^,T) = l2x"(A )W ), (4) 

k 

respectively. We have performed a numerical calculation 
for the integrated spin structure factor (3) and integrated 
susceptibility (4). The result of the integrated spin struc- 
ture factor for the parameter t/J — 2.5 with the temper- 
ature T = 0.3 J (solid line), T = 0.4 J (dashed line), and 
T = 0.5 J (dotted line) at the doping (a) 5 = 0.08, and (b) 
S = 0.15 is plotted in Fig. 2. From Fig. 2, it is shown that 
the integrated spin structure factor is almost temperature 
independent and the shape appears to be particularly 
universal in the underdoped and optimal doped regimes. 
S(uj) is decreased with increasing energies for w < 0.5t, 
and almost constant for ui > 0.5t, which is consistent 
with the experiments 16 and numerical simulations 17 . In 
correspondence with the integrated spin structure fac- 
tor, the result of integrated susceptibility at the doping 
S = 0.15 for the parameter t/J = 2.5 with the tempera- 
ture T = 0.2 J (solid line), T = 0.3 J (dashed line), and 
T = 0.4J (dotted line) is plotted in Fig. 3. For com- 
parison, the function 6iarctan[aiw/T + 03(0; /T) 3 ] with 
&i = 0.23, ai = 2.0, and 03 = 1.4 is also plotted in 
Fig. 3 (dot-dashed line) . Our results show that the inte- 
grated susceptibility is almost constant above ui/T > 1 
and then begin to decrease with decreasing oj/T for 
uj/T < 1. It is quite remarkable that our theoretical 
results of the integrated susceptibility are scaled approx- 
imately as I(uj, T) cx arctan[aiaj/T + a3((jj/T) 3 ], which is 
in very good agreement with the experiments 16 . 

The temperature dependence of the susceptibility con- 
verges to a universal function of uj/T is very significant 
because of its relation to other normal state properties, 
such as the temperature dependences of the spin-lattice 
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relaxation time. The NQR spin-lattice relaxation time 
Ti is expressed as, 

where g is the g factor, /io is the Bohr magneton, and the 
form factors F a (k) = F\\(k)), with F x (k) and 

F|| (k) are for the field applied parallel and perpendicular 
to the C axis, respectively. The form factors have di- 
mension of energy, and magnitude determined by atomic 




FIG. 1. The spinon's second-order self-energy diagram. 
The solid and dashed lines correspond to the holon and spinon 
propagators, respectively. 
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FIG. 2. The integrated dynamical spin structure factor at 
the doping (a) 8 = 0.08 and (b) 8 = 0.15 for t/J = 2.5 with 
the temperature T — 0.3 J (solid line), T = 0.4 J (dashed 
line), and T = 0.5 J (dotted line). 
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FIG. 3. The integrated susceptibility at the doping 
8 = 0.15 for t/J = 2.5 with the temperature T = 0.2J (solid 
line), T = 0.3 J (dashed line), and T = 0.4 J (dotted line). The 
dot-dashed line is the function biarctan[ai<^/T + as(u>/T) 3 ] 
with 6i = 0.23, ai = 2.0, and a 3 = 1.4. 



physics, and k dependence determined by geometry. For 
the comparison with experiments, the form factors F±(k) 
and Fn(k) are chosen as proposed in Ref. 18 . The spin- 
lattice relaxation time T\ in Eq. (5) has been evaluated 
numerically and the results for t/J = 2.5 with the dop- 
ing S — 0.15 for the field applied parallel to C axis (solid 
line) and perpendicular to C axis (dashed line) are plot- 
ted in Fig. 4, where we have chosen units h = Kb = 1. 
From Fig. 4, it is shown that T\ is very weakly dependent 
on T in the optimal doped regime. Some experiments 5 
show that 1/Ti approaches nearly the temperature in- 
dependent at high temperature for I^-ajSr^CuO^ and 
the weakly temperature dependent for YBa2Cu307, in 
the optimal doped regime. Although the simplest t-J 
model can not be regarded as the complete model for 
the quantitative comparison with the copper oxide ma- 
terials, but our results are in qualitative agreement with 
these remarkable experiments 5 . 

In the fcrmion-spin theory, the charge and spin degrees 
of freedom of the physical electron are separated as the 
holon and spinon, respectively. Although both holons 
and spinons contributed to the charge and spin dynam- 
ics, but it has been shown that the scattering of holons 
dominates the charge dynamics 13 , while the present re- 
sults shows that scattering of spinons dominates the spin 
dynamics. The spin dynamics probe local magnetic fluc- 
tuations and are a very detailed and stringent test of 
microscopic theories. Our theoretical results within the 
fcrmion-spin formulism leads to the behaviors similar to 
that seen in the experiments and numerical simulations. 
To our present understanding, the main reasons why the 
present theory is successful in studying the normal-state 
property of the strongly correlated copper oxide materi- 
als are that (1) the electron single occupancy on-site local 
constraint is exactly satisfied during the above analytic 
calculation. Since the anomalous normal-state property 
of the copper oxide materials are caused by the strong 
electron correlation in these systems 1-3 , and can be ef- 
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FIG. 4. The spin-lattice relaxation time at the doping 
8 — 0.15 with t/J — 2.5 for the field applied parallel to C 
axis (solid line) and perpendicular to C axis (dashed line). 
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fectively described by the t-J model 7 , but the strong elec- 
tron correlation in the t-J model manifests itself by the 
electron single occupancy on-site local constraint, which 
means that the electron Hilbert space is severely re- 
stricted due to the strong electron repulsion interaction. 
This is why the crucial requirement is to treat this con- 
straint exactly during the analytic discussions. (2) Since 
the local constraint is satisfied even in the MFA within 
fermion-spin theory, the extra gauge degree of freedom 
occurring in the slave-particle approach does not appear 
here 9,10 , then spinons and holons within the fermion- 
spin theory are by themselves gauge invariant, they are 
real and can be interpreted as physical excitations. In 
this case, there are two relaxation times for spinons and 
holons, respectively, the spinon relaxation time is respon- 
sible to the spin dynamics, while holon relaxation time is 
responsible to the charge dynamics. This important is- 
sue to the copper oxide materials within the charge-spin 
separation has been emphasized by Laughlin 19 . 

In summary, we have studied the spin dynamics of the 
copper oxide materials in the underdoped and optimal 
doped regimes within the framework of the fermion-spin 
theory. Our results show that the strongly correlated 
renormalization effects for the spinon due to the strongly 
spinon-holon interaction are very important for the spin 
dynamics. The integrated dynamical spin structure fac- 
tor, integrated susceptibility, and spin-lattice relaxation 
time are discussed, and the results are qualitative consis- 
tent with the experiments and numerical simulations. 
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